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^ . I. INTRODUCTION 

M\ 

The generalized exponential-cosine-screened Coulomb (GECSC) potential or the generalized 
cosine Yukawa (GCY) potential: 

y (r) = — ( — j exp(— (5r) cos{g5r), (1) 

where A is the strength coupling constant and 5 is the screening parameter, is known to 
describe adequately the effective interaction in many-body enviroment of a variety of fields 
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such as atomic, nuclear, solid-state, plasma physics and quantum field theory [1,2]. It is 
also used in describing the potential between an ionized impurity and an electron in a metal 
[3,4] or a semiconductor [5] and the electron-positron interaction in a positronium atom in 
a solid [6]. The potential in (1) with g = 1 is known as a cosine-screened Coulomb potential. 
The static screened Coulomb (SSC) potential {g = case) is well represented by Yukawa 
form: V{r) = —{aZe'^)exp{—6r)/r which emerges as a special case of the ECSC potential 
in (1) with A = aZe'^, where a = (137.037)""'^ is the fine-structure constant and Z is the 
atomic number, is often used for the description of the energy levels of light to heavy neutral 
atoms [7]. It is known that SSC potential yields reasonable results only for the innermost 
states when Z is large. However, for the outermost and middle atomic states, it gives rather 
poor results. Although the bound state energies for the SSC potential with Z = 1 have been 
studied [7]. 

The Schrodinger equation for such a potential does not admit exact solutions. So various 
approximate methods [8] both numerical and analytical have been developed Hence, the 
bound-state energies of the ECSC potential were first calculated for the Is state using 
numerical [3,8,9] and analytical [10,11] methods and for the s states by a variational method 
[12]. Additionally, the energy eigenvalues of the ECSC potential [13] have been recalculated 
for the Is state with the use of the ground-state logarithmic perturbation theory [14,15] and 
the Pade approximant method. The problem of determining the critical screening parameter 
6c for the s states was also studied [16]. 

It has also been shown that the problem of screened Coulomb potentials can be solved 
to a very high accuracy [17] by using the hypervirial relations [18,19,20] and the Pade ap- 
proximant method. The bound-state energies of the ECSC potential for all eigenstates were 
accurately determined within the framework of the hypervirial Pade scheme [21]. Further, 
the large-N expansion method of Mlodinow and Shatz [22] was also applied to obtain the 
energies of the ground and first excited s— states and the corresponding wave functions. 
Recently, we studied the bound-states of the ECSC potential for all states using the shifted 
large A^— expansion technique [23]. 



In this paper, we investigate the bound-state properties of ECSC potential using a new 
perturbative formahsm [24] which has been claimed to be very powerful for solving the 
Schrodinger equation to obtain the bound-state energies as well as the wave functions in 
Yukawa or SSC potential problem [24,25] in both bound and continuum regions. This novel 
treatment is based on the decomposition of the radial Schrodinger equation into two pieces 
having an exactly solvable part with an additional piece leading to either a closed analytical 
solution or approximate treatment depending on the nature of the perturbed potential. 

The contents of this paper is as follows. In section II we breifly outline the method 
with all necessary formulae to perform the current calculations. In section III we apply 
the approach to the Schrodinger equation with the ECSC potential and present the results 
obtained analytically and numerically for the bound-state energy values. Finally, in section 
IV we give our concluding remarks. 

II. THE METHOD 

For the consideration of spherically symmetric potentials, the corresponding Schrodinger 
equation, in the bound state domain, for the radial wave function reads 

^' ^"(") - V(r) - E (2) 

with 

h^ m + 1 



V(r) 



Voir) + 



+ A\/(r), (3) 



2m r^ 

where AV{r) is a perturbing potential and ipn{r) = Xn{r)un{r) is the full radial wave func- 
tion, in which Xn{r) is the known normalized eigenfunction of the unperturbed Schrodinger 
equation whereas M„(r) is a moderating function corresponding to the perturbing potential. 
Following the prescription of Ref. 24, we may rewrite (2) as 

2m \Xn{r) Un{r) Xn{r)un{r) J 



The logarithmic derivatives of the unperturbed Xn(^) ^^^ perturbed Un{r) wave functions 
are given by 

W'„(0 = --^^ and AW'„=--5=<M, (5) 

which leads to 



^ Xnir) ,,,2,., ^ 



W: r) ^=Wjr) 



Voir) 



h' £{£ + 1) 
2m r^ 



£n, (6) 
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where e„ is the eigenvalue for the exactly solvable potential of interest, and 

f- f ^ + 2 ^"|'|^"|'| ) = AW^ir) - ^A<(r) + 21^„(r)Al^„(r) = AK(r) - Ae„, 

(7) 

in which Ae„ = i^^^^ + E^^ + ■ ■ ■ is the correction term to the energy due to AV{r) and 
En = £n + Asn- If Eq. (7), which is the most significant piece of the present formalism, can 
be solved analytically as in (6), then the whole problem, in Eq. (2) reduces to the following 
form 

[Wn{r) + AWn{r)f - -^(W^„(r) + AW^„(r))' = \/(r) - E„, (8) 

which is a well known treatment within the frame of supersymmetric quantum theory 
(SSQT) [26]. Thus, if the whole spectrum and corresponding eigenfunctions of the un- 
perturbed interaction potential are known, then one can easily calculate the required super- 
potential Wn{r) for any state of interest leading to direct computation of related corrections 
to the unperturbed energy and wave function. 

For the perturbation technique, we can split the given potential in Eq.(2) into two parts. 
The main part corresponds to a shape invariant potential, Eq. (6), for which the super- 
potential is known analytically and the remaining part is treated as a perturbation, Eq. 
(7). Therefore, it is obvious that ECSC potential can be treated using this prescription. In 
this case, the zeroth-order term corresponds to the Coulomb potential while higher-order 
terms consitute the perturbation. However, the perturbation term in its present form cannot 



be solved exactly through Eq. (7). Thus, one should expand the functions related to the 
perturbation in terms of the perturbation parameter A, 

oo oo oo 

AVir;\) = J2\Mr), AW-„(r; A) = ^ A,W-«(r), E«(A) = J] A,i?«, (9) 

j=l i=l j=l 

where i denotes the perturbation order. Substitution of the above expansions into Eq. (7) 
and equating terms with the same power of A on both sides up to O(A^) gives 

2W^{r)Wi'\r) - !-^^f^"^ = V,{r) - E«, (10) 

Vim cir 

Wi^r) + 2W^^{r)W^j^\r) - ^^Wi^) ^ ^^^^^ _ ^^^ ^ 

\j2m dv 

W^W^^) + y^^\rW^;^{r-)\ - l^d^nK'-^ = v,ir) - Ef^. (12) 

-' y 2m cir 

Hence, unlike the other perturbation theories, Eq. (7) and its expansion, Eqs. (10-12), 

give a flexibility for the easy calculations of the perturbative corrections to energy and 

wave functions for the nth state of interest through an appropriately chosen perturbed 

superpotential. 

III. APPLICATION TO THE ECSC POTENTIAL 

Considering the recent interest in various power-law potentials in the literature, we work 

through the article within the frame of low screening parameter. In this case, the ECSC 

potential can be expanded in power series of the screening parameter 5 as [10] 

/ A\ //1\ °° 

V{r) = - - exp(-5r) cos(5r) = - - E '^^i^^)\ (13) 

\r J V ^ / i=o 

where the perturbation coefficients Vi are given by 

V^ = -l,V2 = 0, ^3 = 1/3, V, = -1/6, V, = 1/30, ■ ■ • . (14) 

We now apply this approximation method to the ECSC potential with the angular momen- 
tum barrier 



V(r) = - 



A 



-,—5r 



cos((5r) + 






Voir) + 
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+ AV{r), 



(15) 



where the first piece is the shape invariant zeroth-order which is an exactly solvable piece 
corresponding to the unperturbed Coulomb potential with Vo(r) = —A/r while AV{r) = 
A6 — (A5^/3)r^ + (A5^/6)r^ — {A6^/30)r^ + ■ ■ ■ is the perturbation term. The literature is 
rich with examples of particular solutions for such power-law potentials employed in different 
fields of physics, for recent applications see Refs. [27,28] . At this stage one may wonder why 
the series expansion is truncated at a lower order. This can be understood as follows. It is 
widely appreciated that convergence is not an important or even desirable property for series 
approximations in physical problems. Specifically, a slowly convergent approximation which 
requires many terms to achieve reasonable accuracy is much less valuable than the divergent 
series which gives accurate answers in a few terms. This is clearly the case for the ECSC 
problem [29]. However, it is worthwhile to note that the main contributions come from the 
first three terms. Thereby, the present calculations are performed up to the second-order 
involving only these additional potential terms, which suprisingly provide highly accurate 
results for small screening parameter S. 

A. Ground State Calculations (n = 0) 

In the light of Eq. (6), the zeroth-order calculations leading to exact solutions can be 
carried out readily by setting the ground-state superpotential and the unperturbed exact 
energy as 

mA'^ 



l^„=n r 



' '^' + M^A^, E^ 



'2m r 



2 {e + i)h' 



2h\n + i+iy 



, n = 0,l,2,.... 



(16) 



and from the literature, the corresponding normalized Coulomb bound-state wave function 
[30] 



r(c)^^+i 



Xnir) = K^r'^'exp [-(3r] x Lf+^ [2/3r] 



(17) 



in which A^, 



(C) 



nl 



2mA 



{n+e+i)h^ 



(+1 



1 



;;r^("+2^+i)! 



is a normahzed constant, j3 



mA 



(n+e+i)h^ 



and L^ (x) = E 



m=0\ ''I (n-m)\{m^k)\ra\ 



^"''' '' X™ is an associate Laguarre polynomial function [31]. 



For the calculation of corrections to the zeroth-order energy and wavefunction, one needs 
to consider the expressions leading to the first- and second-order perturbation given by Eqs. 
(10-12). Multiplication of each term in these equations by x^(r), and bearing in mind the 
superpotentials given in Eq. (5), one can obtain the straightforward expressions for the 
first-order correction to the energy and its superpotential: 



E^}'^ 



X^(^) 



Ab'' 



r^ ) rfr, W^'^ (r) 



/2m 1 



Xn(a^) 
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and also for the second-order correction and its superpotential: 






Abr' 



-X 



dx, (18) 



1^(2) 



Ei'^ 
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6 



W^}^' 



dr. 



E^^^ + WJ^'^Hx)-"^^'-' 



6 



-X 



dx 



(19) 



for any state of interest. The above expressions calculate W^^\r) and W,j:^\r) explicitly 
from the energy corrections E^^ and E^^ respectively, which are in turn used to calculate 
the moderating wave function 'U„(r). 

Thus, through the use of Eqs. (18) and (19), after some lengthy and tedious integrals, 
we find the zeeroth order energy shift and their moderating superpotentials as 



E, 
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h\i + if {i + 2) {2i + 3) 
QAw? 



5^ 



E, 
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in which 
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3Am 


Am 


h\^+ir 
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^^(£+l)2(8£2 + 37£ + 43)52 3(2£ + 5) 
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2 (£+1) 
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(21) 



Therefore, the analytical expressions for the lowest energy and full radial wave function of 
an ECSC potential are then given by 



M 



.(1) , 77.(2) 



En=Q/ = El^l^^g + A6 + Ey + Ey -\ , ^„=o,K^) ~ Xn=oA^)'^n=oA^)^ 



(22) 



in which 



Un=oAr) ^ exp 



'2m 



h 



(iVo^^^ (x) + WP (x)) dx 



(23) 



Hence, the explicit form of the full wave function in (22) for the ground state is 



2mA 



(i+lW 



i+i 



Am, 



(£ + i)2\ n\2e + i)\ 



/+^exp(P(r)), 



(24) 



with P{r) = Y^\=iPir^ is a polynomial of fifth order having the following coefficients: 



Pi 



A 



E^^- 



Am 



{^ + i)n- 



u [t -\- Zj 2 1 i-A ^ -fi-A -L r-4 -'- rfi 

'' ^^ " 4(7TT)2^ ' ^^ " 6 ' ^^ " % ' ^^ " 10 ' 

(25) 



in which d = h + rfi((T^\2^ ^-iid other parameters are given in (21). 



B. Excited state calculations (n > 1) 

The calculations procedures lead to a handy recursion relations in the case of ground 
states, but becomes extremely cumbersome in the description of radial excitations when 
nodes of wavefunctions are taken into account, in particular during the higher order calcu- 
lations. Although several attempts have been made to bypass this difficulty and improve 
calculations in dealing with excited states, (cf. e.g. [32], and the references therein) within 
the frame of supersymmetric quantum mechanics. 
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Using Eqs. (5) and (17), the superpotential Wn{r) which is related to the excited states 
can be readily calculated through Eqs. (18) and (19). So the first-order corrections in the 
first excited state (n = 1) are 

QAm? 



3^2^ \ Am \ ^ ' 

Consequently, the use of Wi (r) in the preceeding equation in (19) gives the energy correc- 
tion in the second-order as 

.(2) _ /i'(£ + 2)'(£ + ll)(2£ + 3)(2£ + 5)^4 



h^o (c j-o\^ (c J- 'i\ (Of -L 'X\ (7^^ _L 1 ni ^ -L 91 r, 

'-5\ (27) 



2AA^m? 
h^^ {£ + 2f {£ + 3) {2i + 3) {7f + 101£ + 211) 



72A4m5 

Therefore, the approximated energy value of the ECSC potential corresponding to the first 
excited state is 

En=i,, = Ei'^ +AS + Ei'^ + ^f ) + ■ ■ ■ . (28) 

The related radial wavefunction can be expressed in an analytical form in the light of Eqs 
(18), (19) and (22), if required. The appromation used in this work would not affect con- 
siderably the sensitivity of the calculations. On the other hand, it is found analytically that 
our investigations put forward an interesting hierarchy between W^^' (r) terms of different 
quantum states in the first order after circumventing the nodal difficulties elegantly, 

Wi^\r) = _M!i+l + lH!!: L + n\n + i+l)in + i + 2) \ 
" 3V2m \ Am J ' 

which, for instance, for the second excited state (n = 2) leads to the first-order correction 

^(1) ^ h'ii + 3)'ii + 2)i2i + 23) ^^,^ 
QAm? 

„f.,,) = _M±£)i!:L + ?!!(i±|)(i±i)\ (30) 

^ ' 3^2^ 1 Am \ ^ ' 



Thus, the use of W2 (r) in the preceeding equation (19) gives the energy correction in the 
second-order as 

(2) _ ^' {i + 2) (£ + 3)' (2f + 5) {2e + 45£ + 153) ^ 
^ ~ 24A2m3 

h^° {£ + 2){i + 3f (16£4 + 474£3 + 3879^2 + 12118£ + 12873) ^ 

72lW ^ • ^^^ 

Therefore, the approximated energy eigenvalue of the ECSC potential corresponding to the 

second excited state is 

En=2/ = 4°^ + AS + 4'^ + 4'^ + ■ ■ ■ • (32) 

For the numerical work, some numerical values of the perturbed energies of the Is and 2s 
states, in the atomic units we take h = m = A = 1, for different values of the screening 
parameter 6 in the range < 5 < 0.10 are displayed in Tables 1 and 2, respectively. The 
results are consistent to order 6^ with earlier results obtained by applying different methods 
in Refs. [9,22,23]. Further, we display the results for the energy eigenvalues of 2s, 2p, 3s, 
3j9, and 3d states in Tables 3 and 4. Our results are then compared with accurate energy 
eigenvalues obtained by other authors. Thus, through the comparison of our results with 
those of Refs. [9,10,22,23] for large n and i— values and small screening parameter values 
yields indeed excellent results. 

On the other hand, we take A = ^J2 and b = \plG. Cosequently, we compute the 
binding energies {—En/) of the lowest-lying states. Is to 3d, for various values of 6. Hence, 
the detailed analysis of the results in terms of various domains of parameters A and 6 of 
ECSC potential are displayed in Table 5. For further study of the bound-state energies and 
normalizations with analytical perturbation calculation in Table 6. We consider A = Z = 4, 
8, 16, 24 in order to cover the range of low to high atomic numbers. For low strength of 
A = Z, the energy eigenvalues nobtained are in good agreement with the other methods 
for low values of the screening parameter 6. Obviously, when 6 is small the Coulomb field 
character prevails and the method has been adjusted to that. However, the results become 
gradually worse as A and/or 6 are large. 

10 



IV. CONCLUDING REMARKS 

We have shown that the bound-state energies of the exponential cosine screened Coulomb 
(ECSC) potential for all eigenstates can be accurately determined within the framework of a 
new approximation formalism. Avoiding the disadvantages of the standard non-relativistic 
perturbation theories, the present formulae have the same simple form both for ground and 
excited states and provide, in principle, the calculation of the perturbation corrections up 
to any arbitrary order in analytical or numerical form. 

Additionally, the application of the present technique to ECSC potential is really of great 
interest leading to analytical expressions for both energy eigenvalues and wave functions. 
Comparing various energy levels with different works in the literature we find that this 
treatment is quite reliable and further analytical calculations with this non-perturbative 
scheme would be useful. In particular, the method becomes more reliable as the potential 
strength increases. 
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TABLES 
TABLE L Comparison of bound energy eigennvalues for < 5 < 0.1 for the Is state in atomic 

units. 



6 


1/N [22] 


Dynamical [9] 


Shifted 1/N [23] 


En/ 


0.01 


-0.490 001 


-0.490 001 




-0.490 000 9 


0.02 


-0.480 008 


-0.480 007 8 


-0.480 007 83 


-0.480 007 8 


0.03 


-0.470 026 


-0.470 026 




-0.470 025 9 


0.04 


-0.460 061 


-0.460 060 9 


-0.460 061 01 


-0.460 060 8 


0.05 


-0.450 117 


-0.450 117 4 




-0.450 117 2 


0.06 


-0.440 200 


-0.440 200 4 


-0.440 200 57 


-0.440 200 


0.07 


-0.430 313 






-0.430 313 4 


0.08 


-0.420 461 


-0.420 463 6 


-0.420 463 86 


-0.420 461 7 


0.09 


-0.410 647 






-0.410 648 8 


0.1 


-0.400 875 


-0.400 883 9 


-0.400 884 21 


-0.400 878 5 
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TABLE II. Comparison of bound energy eigennvalues for < 5 < 0.1 for the 2s state in atomic 
units. 



6 


1/iV [22] 


Dynamical [9] 


Shifted l/N [23] 


En,i 


0.01 


-0.115 013 


-0.115 013 5 




-0.115 013 4 


0.02 


-0.105 103 


-0.105 103 6 


-0.105 103 61 


-0.105 103 3 


0.03 


-0.095 334 


-0.095 336 6 




-0.095 334 6 


0.04 


-0.085 755 


-0.085 769 


-0.085 769 59 


-0.085 762 1 


0.05 


-0.076 406 


-0.076 449 7 




-0.076 432 6 


0.06 


-0.067 311 


-0.067 421 7 


-0.067 426 08 


-0.067 390 


0.07 


-0.058 482 






-0.058 680 


0.08 


-0.049 915 


-0.050 392 2 


-0.050 408 25 


-0.050 357 6 


0.09 


-0.041 598 






-0.042 494 5 


0.1 


-0.033 500 


-0.034 967 7 


-0.035 004 67 


-0.035 188 
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TABLE III. Energy eigenvalues as a function of screening parameter S for the 2s and 2p states 

in atomic units. 

State 6 £^[10,10] [10] £[10,11] [10] Pertur.[10] Variational [10] Shifted [23] En/ 

2s 0.10 -0.034 941 -0.034 941 -0.034 425 -0.034 935 -0.035 004 67 -0.035 188 

2p -0.032 469 -0.032 469 -0.032 042 -0.032 470 15 -0.032 673 3 

2s 0.08 -0.050 387 -0.050 387 -0.050 222 -0.050 384 -0.050 408 25 -0.050 357 6 

2p -0.048 997 -0.048 997 -0.048 996 93 -0.048 993 9 

2s 0.06 -0.067 421 -0.067 421 -0.067 385 -0.067 421 -0.067 426 08 -0.067 390 

2p -0.066 778 -0.066 778 -0.066 777 29 -0.066 761 1 

2s 0.04 -0.085 769 -0.085 769 -0.085 767 -0.085 769 -0.085 769 59 -0.085 762 1 

2p -0.085 591 -0.085 591 -0.085 559 13 -0.085 552 

2s 0.02 -0.105 104 -0.105 104 -0.105 104 -0.105 104 -0.105 103 61 -0.105 103 3 

2p -0.105 075 -0.105 075 -0.105 075 -0.105 074 64 -0.105 074 4 
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TABLE IV. Energy eigenvalues as a function of screening parameter S for the 3s, 3p and 3d 

states in atomic units. 

State 6 -E;[10, 10] [10] £;[10, 11] [10] Pertur. [10] Variational [10] Shifted [23] En,e 

-0.005 666 38 -0.007 077 8 
-0.004 492 33 -0.005 405 8 
-0.002 313 56 -0.002 924 
-0.011 685 44 -0.011 952 3 
-0.010 939 85 -0.011 111 7 
-0.009 555 42 -0.009 694 
-0.018 867 16 -0.018 858 6 
-0.018 457 05 -0.018 450 5 
-0.017 682 08 -0.017 691 
-0.036 027 38 -0.036 021 3 
-0.035 967 71 -0.035 964 
-0.035 850 66 -0.035 849 

TABLE V. Energy eigenvalues of the ECSC potential in units oi h = m = 1, A = 2^'^ and 
6 = GA. 



3s 


0.06 


-0.005 461 


-0.005 462 


-0.004 538 


-0.005 454 


3p 




-0.004 471 


-0.004 472 






3d 




-0.002 308 


-0.002 309 






3s 


0.05 


-0.011 576 


-0.011 576 






3p 




-0.010 929 


-0.010 929 


-0.010 538 




3d 




-0.009 555 


-0.009 555 


-0.009 292 




3s 


0.04 


-0.018 823 


-0.018 823 


-0.018 707 


-0.018 822 


3p 




-0.018 453 


-0.018 453 






3d 




-0.017 682 


-0.017 682 






3s 


0.02 


-0.036 025 


-0.036 025 


-0.036 022 


-0.036 025 


3p 




-0.035 968 


-0.035 968 


-0.035 965 




3d 




-0.035 851 


-0.035 851 


-0.035 849 





G State —£^0 o State —Ei q State —Eq i State —Ei i State — £'( 



0,2 



0.002 Is 0.996 000 2s 0.246 000 2 2p 0.246 000 1 3p 0.107 112 3d 0.107 111 4 



0.005 


0.990 000 2 


0.240 003 4 


0.240 002 4 


0.101 125 5 


0.101 116 


0.010 


0.980 001 9 


0.230 026 9 


0.230 019 3 


0.091 221 7 


0.091 147 5 


0.020 


0.960 015 6 


0.210 206 6 


0.210 148 9 


0.071 928 1 


0.071 361 7 


0.025 


0.950 030 2 


0.200 395 3 


0.200 285 7 


0.062 648 5 


0.061 566 5 


0.050 


0.900 234 4 


0.152 865 2 


0.152 099 1 


0.022 223 5 


0.014 137 4 



TABLE VI. Energy eigenvalues of the ECSC potential for all states in units of 7i = 2m = 1, 
and S = 0.2 fm^^. 

A i n —En/ A i n —En/ 
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3.207 029 








14.403 752 


1 





2.433 587 








60.801 938 


1 





12.818 287 








139.201 31 


1 





31.212 563 


2 





11.249 961 



16 60.801 938 24 



24 






1 


12.825 303 





2 


4.023 139 


1 


1 


4.009 505 





1 


31.217 455 





2 


11.279 786 


1 


1 


11.269 899 


1 


2 


4.412 177 


2 


1 


4.380 887 


2 


2 


1.411 568 
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